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Abstract The optical spectral energy distributions of two tidal disruption flares identified by van Velzen et 
al. (2011) in archival SDSS data, are found to be well-fit by a thin-accretion-disk model. Furthermore, the in- 
ferred Supermassive Black Hole mass values agree well with the SMBH masses estimated from the host galaxy 
properties. Integrating the model SEDs to include shorter wavelength contributions provides an estimate of the 
bolometric luminosities of the accretion disks. The resultant bolometric luminosities are well in excess of the 
minimum required for accelerating UHECR protons. In combination with the recent observational estimate of 
the TDF rate (van Velzen and Farrar, these Proceedings), the results presented here strengthen the case that 
transient jets formed in tidal disruption events may be responsible for accelerating all or most UHECRs. 



1 Introduction 

It was proposed in ref. ^ (FG09 below) that all or most 
UHECRs may be accelerated in transient jets produced by 
stellar tidal disruption events or extremely powerful disk 
instabilities in AGNs. This proposal was motivated in part 
by the shortcomings of GRBs and powerful AGNs, the 
primary contenders for UHECR acceleration, and in part 
by the association reported by the Auger collaboration (O 
between UHECRs and galaxies in the Veron-Cetty Veron 
(VCV) list of possible AGNs(3 ). Of Auger's first 27 events 
above ~ 55 EeV iS, 19 of the 22 events with \b\ > 10° 
are within 3.1° of a VCV AGN candidate (excluding the 
Galactic plane where galaxy catalogs are not complete). 

FG09 pointed out that such a correlation with AGNs 
would be surprising, unless UHECR production occurs in 
a transient state in which the luminosities of the jet, and the 
accretion disk powering it, are much higher than found in 
all but a few of the most powerful steady-state AGNs. The 
argument goes as follows. To correlate within 3.1° implies 
that most UHECRs are protons, because high-Z heavy pri- 
maries would be deflected much more than a few degrees 
in the Galactic magnetic field, and intermediate-Z UHE- 
CRs would photodisintegrate enroute from the source. A 
proton can only be accelerated to energies such that its Lar- 
mor radius is smaller than the size of the accelerator, plac- 
ing a lower bound on BxR, the magnetic field times source 
size. This implies that the (isotropic equivalent) total power 
in the required magnetic field (Poynting luminosity) is of 
order 

L ~ IcI^B^R^ > IQ'^^r^ {E2QlZf erg/s. (1) 
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If the energy in the magnetic field, protons and electrons 
is in equipartition, and the energy in electrons is emitted 
through synchrotron cooling in the time it takes the shock 
to pass through the magnetic cloud, a comparable lumi- 
nosity is also emitted by electrons with jet-frame Lorentz 
factors lO-' <ye < 10^ 5 placing a limit on the luminos- 
ity of the jet. In order to maintain such a jet, the accretion 
power should be at least as high, and the bolometric lu- 
minosity of the AGN should satisfy Eq. [T] GRBs easily 



satisfy this limit, but only the most powerful AGNs satisfy 
itforZ = 1. 

Neither GRBs, powerful AGNs, nor arbitrary sources 
embedded in the local large scale structure (distribution of 
galaxies as determined from the 2MASS redshift survey) 
would give rise to a UHECR-AGN correlation on small 
angular scales at the level reported(5 ). Furthermore, FG09 
showed from the observed GRB rate and luminosity distri- 
bution, that GRBs can only account for the total UHECR 
flux if the power they emit in UHECRs is much larger than 
the power emitted in gammas, which would be surprising 
theoretically. FG09 proposed instead that UHECRs are ac- 
celerated in jets with a lifetime of order months, produced 
by TDEs or exceptional accretion disk instabilities. An as- 
sociation with AGNs would arise if the presence of a thin 
accretion disk enhanced the probability of a star being cap- 
tured, or enhanced the probability of the accretion process 
producing a sufficiently powerful jet. FG09 showed that 
given the predicted rate of TDEs ^ and a 1% efficiency, 
TDFs could easily satisfy both the peak luminoisty require- 
ment, Eq.[T] and the total UHECR energy injection require- 
ment. 

In the following, the SEDs of the two TDFs found in 
SDSS Stripe 82 archival data by van Velzen et al. (2012) 
([7]l (vVFl 1 below), are used to infer Lboi for these TDFs, 
to learn if TDFs do indeed satisfy the peak luminoisty re- 
quirement, Eq.[T] In another contribution to this workshop, 
S. van Velzen and the author estimate the TDF rate us- 
ing the vVFl 1 events. It is found that TDEs satisfy both 
requirements for UHECR production: peak luminosity of 
individual events and adequate total energy injection. 

Before turning to the main task, we note these addi- 
tional, relevant pieces of information: 

1. Zaw, Farrar and Greene(^) followed up on the VCV 
galaxies correlating with UHECRs in the first Auger 
data release. They found that some of the VCV galax- 
ies were not AGNs, some needed additional observa- 
tions to classify, and of the ones that are AGNs, most 
were too weak to satisfy Eq.[T] 

2. Terrano, Zaw and Farrar(9) performed Chandra obser- 
vations on the correlating VCV galaxies whose status 
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as having an active nucleus remained uncertain, and on 
one VCV AGN whose bolometric luminosity could not 
be determined from existing observations(8 ); they also 
classified all but one of the VCV galaxies correlating 
with UHECRs in the second Auger data release (fTOt 
and found Lboi for all but two cases. 
Combining both Auger published UHECR datasets, there 
are 57 events with \b\ > 10°. Ref. (9) found that 30- 
50% of these events correlate, within the 3.1° prescribed 
by Auger, with genuine AGNs with z < 0.018. Of 
these correlating AGNs, two satisfy Eq. [T[ two do not 
have measured Lboi, and the rest do not satisfy Eq. [T] 
Since about 45% of the sources of protons above 55 
EeV should have z < 0.018 due to the GZK effect, the 
observations are consistent with somewhere between 
roughly half and all of the UHECRs being protons cor- 
relating with weak AGNs. 

3. A cluster of 4 UHECRs in the combined AGASA and 
HiRes data (11) (the "Ursa Major" cluster) is consis- 
tent with being protons produced in a transient source 
STli . The small angular dispersion of the cluster im- 
plies the UHECRs are protons and the absence of a 
noteworthy source candidate in the field (TSl argues 
that the source is transient. Furthermore, the UHECR 
energies show the peaked distribution characteristic of 
a transient source and not the power-law distribution of 
a continuous source (fT4l[T5T l. 

4. Auger has measured the distribution of shower max- 
ima, Xmax, in hybrid events up to about 30 EeV. Taking 
the predictions of available shower simulations at face 
value, the observations suggest a shift to heavier com- 
position with increasing energy. However the simula- 
tions fail to describe the observed muon content of the 
ground showers correctly (16), so it is not yet possible 
to draw conclusions from the X^ax distribution; data is 
also needed at higher energy, because it is the compo- 
sition of the UHECR dataset which is in question. 



2 Fitting the TDF SEDs 

The reader is referred to vVFl 1 for details of the analysis 
of SDSS Stripe 82 observations which lead to the iden- 
tification of two tidal disruption events called TDEl and 
TDE2, with redshifts z - 0.136 and 0.251 respectively. 
Their estimated BH masses are (0.6-2) and (2-10) xlO^^°^Mo, 
where the +0.3 in the exponent reflects the scatter in the 
bulge-BH mass relationship, and the range in prefactors 
comes from the uncertainty in the bulge masses of the host 
galaxies. Figs.[T]and|2]show the SEDs for the two TDEs in 
the rest frame of their respective BH; the relatively small 
error bars result from combining observations on different 
nights by rescaling their mean to the initial observation, 
there being no significant color evolution over the 3 months 
of observations. 

For such massive BHs as TDEl, 2, the super-Eddington 
phase is very short, so the observations most likely occured 
in the accretion-disk dominated phase. Lodato and collab- 
orators SVT\ [TSJl and Strubbe and Quataert ( fT9] l discuss the 
light-curves of TDEs and derive the accretion-disk equa- 
tions in the thin-disk approximation, for an assumed M. 
They find, as one would expect, that the emission is a su- 
perposition of black-body annuli, with temperature drop- 
ping with radius. Because the SDSS observations are in the 
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Figure 1: Thin accretion disk fit to the SED of TDEl; data 
from vVFl 1. 
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Figure 2: Thin accretion disk fit to the SED of TDE2; data 
from vVFl 1. 



optical, the observed SEDs are dominated by large radii in 
the accretion disk. 

To restrict the number of free parameters in this ini- 
tial investigation, the disrupted star was assumed to have 
the mass and radius of the Sun and default values of ac- 
cretion efficiency, outflow fraction, etc were adopted from 
(IT9] l; the ratios of pericenter to tidal disruption radii were 
set to their most probable values for the given Mbh^- The 
most important parameters to the fit are the observation 
time, Mbh, and outer radius of the accretion disk; these are 
allowed to vary. 

The quality of both fits is excellent, as can be seen in 
Figs, [l] and |2] The best fit values ai-e Mbh = 1-9 10^ Mg 
and 5.6 10^ Mq, for TDEl, 2 respectively, in good agree- 
ment with the values inferred from the host galaxies. The 
outer radii of the accretion disks are inferred to be 155 
and 205 in units of /?sch- The inferred observation epoch 
is about 6 months for each, consistent with the 9-months 
off, 3-months on obsevation schedule for Stripe 82 and the 
fact that the flares were both first seen at the beginning of 
an observing season. Integrating over all frequencies gives 
bolometric luminosities at the time of observation of 5 and 
10 X 10"*^ erg s"', respectively. 



3 Conclusions 

We have found that a simple model of thin-disk accretion 
in a tidal disruption event, based on formulae derived in 
ref. ( |T9] l, gives a good fit to the observed optical spectral 
energy distributions of both SDSS tidal disruption events. 
If this model is valid, only a small fraction of the total lu- 
minosity is emitted in optical wavelengths, and Lboi easily 
satisfies the criterion for accelerating protons to 100 EeV 
and above. 
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The model does not take into account the temporal evo- 
lution of the accretion disk, and numerical simulations will 
be needed to obtain a more realistic description. Of partic- 
ular interest is how plausible it is for the accretion disks 
to extend to the large radii needed to produce the observed 
optical emission within this model. Could such a large ra- 
dius be a hint of a weak, pre-existing AGN? (The vVFl 1 
analysis excludes AGNs based on locus in a color-color 
plot and on variability, but of course a sufficiently tenuous 
accretion disk could go undetected. 

Much more work is needed in modeling the observed 
SED. Some of the next steps, in addition to developing 
more detailed models of the emission, will be to 1) ex- 
ploit GALEX and Catalina Real-time Transient Survey ob- 
servations of these two TDEs, which have not been used 
here, 2) explore more extensively the space of model pa- 
rameters in the simple thin-disk model, and 3) place up- 
per and lower limits on Lboi. Another interesting question 
is whether a unified model can be developed which de- 
scribes not only these two SDSS events, but also the recent 
Swift events - presumably TDF's viewed from an angle 
such that jet emission dominates emission from the accre- 
tion disk( 20112 lb. 

A clearer picture should emerge from UHECR shower 
observations over the next few years, as to the nuclear com- 
position of UHECRs. Then, we will know better how strin- 
gent the luminosity requirement ([T]) on the sources really 
is. If a substantial fraction of UHECRs are protons or low- 
Z nuclei, bolometric luminosities in excess of 10"*^ erg s"^ 
will be required and at the same time a rare-source scenario 
will be ruled out, unless extragalactic magnetic fields pro- 
duce much larger deflections than the Galactic field does; 
at least a portion of UHECRs must in this case be acceler- 
ated in transient sources, and Tidal Disruption Flares will 
be a strong contender. As more cases like the Ursa Major 
UHECR cluster with multiple events from a single source 
are found, whether or not the source is bursting can be de- 
termined by the characteristic peaked spectral shape. 
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